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We have developed a theory of the negative magnetoresistance due to the weak localization in 
uniaxial wurtzite-type crystals, in which the spin splitting of the conduction band is linear in the 
wave vector, unlike the cubic A3B5 crystals. Unlike earlier theories, we take into account the 
correlation between the electron motion in spin and co-ordinate spaces. It is shown that as a result 
of this correlation the magnetoresistance depends on the orientation of the magnetic field with 
respect to the main axis of the crystal even if the effective mass is isotropic. The new theory allows 
to accurately determine the value of the spin splitting constant in uniaxial crystals. 
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The effect of the negative magnetoresistance observed 
in highly doped semiconductors and metals in weak mag- 
netic fields is known to be caused by the weak localiza- 
tion. The weak localization, which is the result of the 
constructive interference of two electron waves traveling 
along a closed path in opposite directions and scattering 
by the same impurities, leads to an additional contribu- 
tion to the resistance. In a magnetic held the two waves 
acquire a phase difference ^/^o, where $ is the magnetic 
flux through the electron trajectory and <&o = fic/2e is 
the flux quantum. As a result, the interference condi- 
tions are violated and the additional resistance decreases, 
which manifests itself as a negative magnetoresistance, 
i.e. increase of the conductivity JO. 

If one takes into account the electron spin, the sin- 
glet state of the two waves with the total momentum 
J = gives a negative contribution into the resistance 
(antilocalization), while the triplet states with J = 1 
give a positive contribution S. In certain situations 
the triplet contribution is suppressed by the spin relax- 
ation, while the singlet contribution is not. Then a weak 
magnetic field first reduces the negative contribution of 
the singlet state, which results in a positive magnetore- 
sistance. The negative magnetoresistance takes over in 
larger fields, which suppress the triplet states stronger 
than the spin relaxation does. According to |^,|| this 
change of the sign of magnetoresistance in bulk crystals 
can be observed for Elliott-Yafet and Dyakonov-Perel 
mechanisms of the spin relaxation. On the other hand, 
the spin relaxation dominated by the scattering on para- 
magnetic impurities suppresses the singlet contribution, 
and the magnetoresistance is negative in the whole range 
of magnetic fields where the weak localization is impor- 
tant. Therefore, the type of the magnetoresistance is 
determined by the dominant spin relaxation mechanism. 
In the non-centrosymmetric semiconductors (with the ex- 



ception of narrow-gap and gapless semiconductors) this 
is the Dyakonov-Perel mechanism, which is caused by 
the spin splitting of the conduction band . 

In the uniaxial wurtzite-type crystals, unlike the cubic 
crystals A3B5, the spin splitting of the conduction band is 
linear in the wave vector fc, as described by the following 
Hamiltonian 101: 



H = a[ak] z 



(1) 



where z is the main axis Cq of the crystal. It was shown 
0-^] that in this case the theory of the weak localization 
must take into account the correlation between the elec- 
tron motion in co-ordinate and spin spaces, similarly to 
the 2D structures. 

The weak localization contribution to the conductiv- 
ity can be expressed through the Cooperon C as follows 

&# 



Act = -2e 2 Di/ T ( 



o 2 E/ 



C a (3f3a{q) 



d 3 q 
(2tt) 3 



(2) 



\v 2 ti is the diffusion coefficient, t\ is the 



Here D = 

transport time, tq is the momentum lifetime, v is the 
electron velocity, g^ax = V > 

v a = (2m) 3 / 2 y/eF I (2ir) 2 is the density of states at the 
Fermi level tp, and m is the effective mass. The summa- 
tion in Eq. (^J) is carried over spin indices a, (3 = ±1/2.Q 
For simplicity we assume an isotropic electron spectrum. 
Then the matrix equation for the Cooperon can be writ- 
ten similarly to Refs. 0-0J: 



* Throughout this paper we take % = 1 except in the final 
expressions 
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Here, according to Eq. (|l|), 

£l x (g) = f2sini?sin<p, 

£ly{g) — — Q sini? cos ip, Q = akF, v g 
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and t v is the phase-breaking time. In Eq. (||) the Pauli 
matrices a act on the hrst pair of spin indices, while 
the matrices p — on the second pair. The angles $ 
and ip are that of the vector v in the coordinate sys- 
tem with x || Cq, and the integration is over solid angle 
4ir: dO g = sint? d§d<p. Vt ?, is the scattering matrix ele- 
ment, which includes the density of the scatterers. Then, 
for the scattering probability we have 



W(k, k') = 2irv 



If the scattering probability depends only on the scat- 
tering angle 9, then the equation 



dO a 
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has the following harmonics as its eigenfunctions: 

C2J = C n>ro Y£0?V)- (5) 

These formulas are which is written in the coordinate 
system with z' \\ is', and and <fr' are the polar and 
axial angles of k in this system. The following property 
is well-known for the spherical harmonics: 

/ ^^ Y ™(0 9 )=Yr«(Ofc) W(0) cosnO sine d6 
= Y^(O k ) (---). (6) 

Here 
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W{9) (1- cos n0) sin 6d6, (n ^ 0), 



ir/2 



W{6) sine dO. 



(7) 



it follows from that the eigenfunctions (||) have the 
corresponding eigenvalues X n , m — 1 — to/t„, (n ^ 0), 
and Aq = 1. Therefore, the only large component of 



the Cooperon is the scalar Co,o = Co- Since the right- 
hand side of Eq. contains linear in g terms, it is also 
necessary to retain the harmonics Cj m . These can be 
found from (|j) if we substitute Cg g = Ci ltn in the right- 
hand side and g — Co + Ci, m in the left-hand side, 
and keep only terms linear in g. Then, using Eqs. (|| 0) 
we can express Ci through Cq: 



(wfe9)To + (a + p)Qt Cq. (8) 
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We now substitute C = Co + Ci to the right-hand side 
of Eq. (H), express Ci through Co using Eq. (||), and 
integrate both sides over dOk and dOk> to cancel out 
all components C n ,m except Co- We then arrive to the 
following equation for the matrix S = 2itvqTq Cq: 



HS = 1, 



(9) 



where 



H = Dq 2 + {d + p) 2 Q 2 n + 

(2Dr 1 ) 1 / 2 n [{a x + p x )q y - {a v + p y )q x ) + — . (10) 

t v 

The Green's function for the equation (|l?j) can be writ- 
ten as [] 



(ii) 



where \I/^g(g) and E r (q) and the eigenfunctions and 
eigenvalues of the operator 7i: 



H^ r (q)=E r ^ r (q). 



(12) 



The basis of eigenfunctions \& r can be converted by the 
appropriate choice of their linear combinations to a basis 
of singlet and triplet states. The Hamiltonian (|lFJ) is then 
split into two matrices, since it has no non-zero matrix 
element between singlet and triplet functions. For the 
singlet state with J = and antisymmetric wave function 
the Hamiltonian is simple: 



rl (q) = Dq 2 



1 



(13) 



while for the triplet state with J = 1 and symmetric wave 
functions ^>} n (m = 0, ±1) it is convenient to introduce 



^Note that the formula ( p"lj ) is valid only when the operator 
7i is hermitian. It is, therefore, inapplicable, for example, 
when the Zeeman splitting in the magnetic field is taken into 
account. 
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the momentum operator J = (a + p)/2 and write the 
Hamiltonian in this form: 

Hi(q) = Dq 2 + 2 (2 - J 2 ) n 2 n 

+2i(Dt 1 ) 1 / 2 Q (J+q_ - J_q+) + — , (14) 



where 

q± = q x ± iq y and J± = —j= (J x ± J y ) ■ 

Since the eigenfunctions f r are orthonormal, it follows 
from Eq. ((jlj) that 



(15) 



Due to the anisotropy of the spin splitting ([!]) the Hamil- 
tonian Tii in a magnetic field depends on the field direc- 
tion. In a magnetic field B \\ z we can use the commuta- 
tion relation 



| [q x q v ] = k+q-} = jj, 



where 
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to introduce new operators 
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(16) 
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The Hamiltonian (|| |J) can be then written as 



Hofe) = £<? 2 + <W} + — , 

H\{q z ) = Dq 2 z + S{aa^} + — + 2 (2 - J 2 ) Q 2 Tl 

+2i(D Tl ) 1/2 n (j+at - J_a) , (19) 

In the basis of the eigenfunctions of the operator {aa< } = 
^(acr + a* a) these operators have the following non-zero 
matrix elements 



(n — 1| a\ri) = (n| a) |n — 1) = 
(n| {aa^} |n) = n + — . 



(20) 



According to Eqs. <M), ( |15|) the weak localization cor- 
rection to the conductivity can be written as 

A<t( j B)=-0 / ^ dq z S{q z ), (21) 



where 



s (v*) = E - 



E 



n=0 



(22) 



E™(q z ) are the eigenvalues of Tto and Wi and 
1/8t\. For B || z these eigenvalues can be calculated an- 
alytically, similarly to Ref. ]7j]. Since we are only inter- 
ested in the change of the conductivity with the magnetic 
field, we subtract Aa(B = 0) from Eq. ( [H|) to obtain the 
magnetoconductivity Scr(B) = Acr(B) — Acr(0). In 5a, 
following Rcf. [[To| , we can extend both the integration 
over q z and summation over n to infinity: 



5a(B) 
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F{x)-- 
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(23) 



1 

ao ' ai(ao + 6) — 26 

3a 2 +2a„6-l-2(2n+l)& 



El 

n=l 
2 



Ofi-l a n-l + b 



1 



1 



26[(2n+ l)a r . 



a n (x) 



Bu- 



ck 



4eDr ir . 



~B~^' 

Hso = 



+ x* 



Hso 
B 

n 2 Tl ch 

2eD 



so 



B 

3a 2 m 2 c 
2eh 3 ' 



x = q z l, (24) 



*»(*) = £ 



n + 1 + z + yTi + Z 
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and / = (ch/eB) 1 / 2 is a magnetic length. 

For Biz we can choose axis x \\ B we can introduce 
o and similarly to Eq. (|l7j): 

£>(«2 + «») - Haai}, 2D 1 / 2 q y = ziS 1 / V _ a ). (25) 
Then the Hamiltonian (141) becomes 



1 



Wife) = Dg£ + £{aa T } + — + z(2(5r 1 ) 1 / 2 r!(a T - a) J. 



2{2DT 1 ) 1/2 VLq x J l 



(26) 



In this case the eigenvalues of the Hamiltonian can be 
computed only numerically, restricting the matrix size to 
some large, but finite, value. It is convenient to express 
the sum of the reciprocal eigenvalues E™(q x ) through 
the minors \Da(q x )\ of the diagonal elements Hi M and 
the determinant \D\ of Hi as follows: 



E E 
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\Dii(q x )\ 
rlj E™( qx ) - ^ \ \D\ 

X>1„(&)- 



(27) 
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FIG. 1. Theoretical magnetoconductivity Sa(B) for B _L z 
(curve 1) and B \\ z (curve 2) for H v = 3Gs and 
//so = 445Gs. The curve 3 shows the magnetoconductiv- 
ity as obtained from the theory of Ref. Q Eq. (^8|) for the 
same values of the parameters. 
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FIG. 2. Comparision of the measured in Ref. |H| (circles) 
and theoretical magnetoconductivity. Solid line show the re- 
sult of our theory for H v = 3 Gs and //so = 445 Gs; dashed 
line is the best fit fit obtained from the theory of Ref. 
(H v = 1.7 Gs and Hso = 141 Gs). 



The first sum converges rapidly for i — ► oo. The summa- 
tion in the second sum can also be extended to infinity if 
we subtract its value at B = similarly to Eq. (|23|). 

If we omit the linear in f2 and q (or a, a') terms in 
Eqs. (£§) and (||), the first terms in Eqs. (§|) and (§7|) 
also vanish. In this approximation the magnetoconduc- 
tivity 5a does not depend on the orientation of the mag- 
netic field and is given by this formula: 



Sa{B) 



2/3 
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(28) 



This formula differs from the one derived in Ref. |3|] 
and used in, for example, Ref. Jll[] in the value of 
H SO - instead of H S o = 1/5t sxx with Ur sxx = l/T syy = 
1/2t szz = 2r» 2 T! Hi, the Rcfs. |,0 used 1/t sxx = 

In Fig. |l| the magnetoconductivity 5a(B) as obtained 
using Eqs. (|l9| - ||), (|e|), and (§?]) is shown for two 
orientation of the magnetic field B \\ z and B _L z for the 
same values of H v = 3Gs and Hso — 445Gs. The dashed 
line shows 5a from Eq. ( p8| ) for the same values of the 
parameters. One can see that 5a significantly depends on 
the field direction and differs strongly from the result of 
Eq. ( |28| ) , which does not take into account the correlation 
of the electron motion in spin and co-ordinate spaces. 

The results of magnetoconductivity measurements for 
CdSe were published in Ref. [O. The samples used had 
electron density n = 8.6 x lO 1 ^ cm~ 3 and kpl = 2.62, 
which corresponds to D = ^ — k F l = 7.84 for 

* 6 m c sec 



m = 0.13 mo, mo being the free electron mass. In Fig. Q 
we reproduce the experimental results of Ref. |nj for 5a 
at T = 0.06 K and compare it with the results of our 
theory. In Ref. Jll[ the direction of the magnetic field 
is not specified; it is noted that at T = 1.5 K 5a(B) 
does not depend on the field orientation (at this temper- 
ature H v > Hso and the dependence 5a(B) does not 
have a minimum at small B, and the role of the linear 
terms in Eqs. (jljt), ( pfj| ) becomes non-essential). The 
solid line in FigTg shows the fit by the theoretical mag- 
netoconductivity (Eqs. (|| - H)). for B || z. The fit- 
ting was done by weighted explicit orthogonal distance 
regression using the software package ODRPACK |l|. 
The weights were selected to increase the importance of 
the low-field (B < IkGs) part of the magnetoconductiv- 
ity curve. The parameters of the best fit are H v = 3 Gs 
and Hso = 445 Gs (these are the same values we used 
in Fig. |l|). The dashed line shows the best fit obtained 



from Eq. (|28[), the parameters of this fit are H v = 1.7 Gs 
and Hso = 141 Gs. One can see that the agreement with 
experiment is not as good as that of the solid curve. 

The above value of Hso — 445 Gs corresponds to 
a = 3.9 x 10~ 2 eV • A. The standard deviation of a as 
computed by the fitting procedure is 3%. In Ref. |ll[] 
authors have obtained a = 2.5 x 10" 2 eV-A (they 
have used the notation A = 2a for the spin splitting 
coefficient). Taking into account the error in the def- 
inition of t sxx pointed above this result should be di- 
vided by V2- In Ref. |l3) the constant A was measured 
for Cdi_ x Mn x Se using electric-dipole spin resonance and 
found to be a = A/2 sa (3.25±1.2) x 10~ 2 eV • A (extrap- 
olated to x = 0). However, for more accurate comparison 



with our theory we must know the precise orientation of 
the magnetic field. We should also note that the weak lo- 
calization theory includes only the corrections to the con- 
ductivity which are of the first order in h/riep — 2/k-pl. 
For the sample used in Ref. jll[] 2/fcpZ = 0.7, which means 
that the theory should not give very accurate results. 

In conclusion, we have developed a theory of the 
weak localization and the negative magnetoresistance in 
wurtzite-type semiconductors. It is demonstrated that 
if the spin splitting of the conduction band is linear in 
the wave vector it is necessary to take into account the 
correlation between the electron motion in co-ordinate 
and spin spaces. This correlation and the anisotropy of 
the spin splitting result in the mangetoresistance being 
dependent on the magnetic field direction even for an 
isotropic electron effective mass. 
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